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Abstract 

We describe the classical two dimensinal nonlinear dynamics of cold atoms 
in far-off-resonant donut beams. We show that there chaotic dynamics exists 
for charge greater than unity , when the intensity of the beam is periodically 
modulated. The two dimensional distributions of atoms in (x, y) plane for 
charge two are simulated. We show that the atoms will acumulate on several 
ring regions when the system enters to regime of global chaos. 
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I. INTRODUCTION 



The Lagurre-Gaussian beam carries orbital angular momentum associated with helical 
surface of constant phase has recently been the subject of considerable theoretical and exper- 
imental study [1-7]. Many of these studies concern the transfer of orbital angular momentum 
of the beam to small particals even atoms . With the developement of laser cooling and 
trapping of neutral atoms, various schemes to slow the atomic motion and increase the 
atomic intensity have been proposed and some demonstarted. In a far-off-resonant laser 
beam it is possible to observe the spatial variation of the optical dipole force on slow atoms. 
In this case the atom experiences an effective mechanical potential propotional to the inten- 
sity of the beam. Recently the Gaussian-Laguerre modes were used as as optical potential 
to trap atoms |J. 

The degree of nonlinearity of the potential is characterised by the 'charge' of the donut. 
The charge is an integer characterising the phase singularity of the beam on axis and which 
also determines the variation of the intensity off axia in the transverse plane. In the case of 
unit charge the intensity close to the axis varies quadratically with the radical distance in 
the transverse plane. For a charge greated than inity the variation of intensity is a quartic, 
or higher, power of the radical distance and thus the resulting motion near the beam axis is 
nonlinear. In this case the frequancy of bound oscillatory motion near the axis depends on 
the energy of the atom. As the energy increases, the frequency decreases, eventually falling 
to zero quite quickly as the unstable fixed point at the potential maximum is reached. It 
is well known that if such a nonlinear system is driven by a periodic modulation of the 



potential, chaotic dynamics can result ||11|| . In this paper we will consider the motion of cold 
atoms in intensity modulated Gaussian-Lagurre (donut) beam and discuss the classical two 
dimensional chaotic dynamics for charge 1 = 2. 
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II. LAGUERRE-GUASSIAN BEAM AND DONUT BEAM 



The expression of a linearily poralized L-G beam in cylindral coordinator(r, 0, z) is 



l + z 2 /Z R 2K w(z)' q y w 2 (z) J 
exp(-r 2 /w 2 (»)exp[ ]exp[i/0] x 

exp[i(2g + I + VjtarT 1 ^/ zr)] exp(ikz), (1) 

where z R = 7£1 f^-, w 2 (z) = ^3^; I = n — m and g = min(n, m) for mode. 

We only concentrate on the donut beam, because it is the easiest and the most applicable 
way to form the optical trap, in which case it has q = and therefore LqQ = 0. In reality 
(1) can be greatly simplified because in general z <^ Zr and w(z) ~ wo. 

The new expression of electric field for donut beam is 

r\/2 

E(r) = e x C(——y 11 exp(-r 2 /wo) exp[i/0] exp(ikz), (2) 

Wo 

where I is the topological charge of singularity for donut beam and the sign can be positive 
or negative. Each photon in the donut beam carrys L z = lh orbital angular momentum for 
linearly polarized donut beam. The intensity distribution of the LG q=0 beam is given by 

2l»l+i r 2|i| 

Hr) = W exp[-2r 2 / Wo 2 ], (3) 

7l\l\\Wo 

where W is the power of laser. 



III. CHAOTIC DYNAMICS OF COLD ATOMS IN FAR-OFF-RESONANT 

DONUT BEAM 

For cold atoms in far-off-resonant donut beam, the spontaneous emission can be ignored. 
The effective optical potential for two-level atom || has the form 

U(r) = ^ln(l+p), (4) 
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where A is the detuning and p = ^^2/4 is a saturation parameter, where Q is Rabi 
frequency, for for-off-resonant donut beam p <C 1 and thus 

Taking into account the special variation of the Rabi frequency, the classical Hamiltonian 
in (x, y) plane for the system is 

#o = 4^ + k{W + v 2 ) 1 (6) 

where K(l) = ^ I ; ,^ 2i and f2o = 2ww 2 f i ^ s * s ^ ne saturation intensity, for Rb atom 
I s = 2mW/cm 2 . 



When the Rabi frequency is modulated and f2 is substitued by f2oy 1 + ecos(cjt),, the 
dynamics of the atom is chaotic . For the different charge /, the potential is different. The 
larger the value of /, the wider the bottom of the potental. For I — 1, at the bottom of this 
potental the motion of atom is harmonic and no chaotic dynamics for periodic modulated 
potential can arise. Only if the atom has a larger kinetic energy so that the nonlinear parts 
of the potential are expected, does the dynamics became chaotic. For I > 1, the dynamics 
is nonlinear over the whole potential range and because the bottom is wider, it is easier to 
control the amplitude of modulation. 

For simplicity we discuss I = 2 case and the practical parameters for our numerical 
examples are that for Rb atoms, the linewidth T/2n = 6MHz, mass M = 85m p , the beam 
waist of laser w = 140/xm, laser power W = 600mW/cm 2 and the detuning A/2n = QGHz. 
We define dimensionless parameters (x,y) = (x,y) = (x/w,y/w), (p~ x ,Py) = (Px/PD,P y / Pd), 
and H = H/2Eq and t = t/jfp^, where E D = is the Doppler limit energy and Pd is the 
momentum respectively. Omitting the tildes and defining £ = 2 ^e d > ^ ne Hamiltonian can 
be rewritten for charge I = 2 as 

2 2 

H(t) = + £(/)(x 2 + y 2 ) V 2( * W) (1 + ecoscut), (7) 

where £ ~ 0.887. Using Hamiltons equations we find that the motion in the transverse plane 
without modulation (e = is described by the equations, 



p' x = — 4£:rr 2 (l — r 2 )e 



(8) 



p y = -4£yr 2 (l-r 2 )e- 2r \ (9) 
x=p x , (10) 

y = Pv, (ii) 

where r 2 = x 2 + y 2 . Clearly there are twofixed points, one stable fixed point on axis (r = 0), 
and one unstable fixed point at the intensity maximum (r=l). 

The chioce of modulation frequency lj depends on the frequency of unperturbed periodic 
motion. For simplicity we assume y = and p y = 0, so the expression for H is simplified as 
one dimensional Hamiltanian. The motion period for unperturbed Hamiltanian H is 

T = I dx = 2 f " — (12) 

J dH /dp x J-xm h(H - £x 2 exp (— 2x 2 ) 



where xm is determined by H = ^x 2 M exp(— 2x 2 M ). Therefore 

7T 

^ = r^{2[^o-^I f exp(-2x 2 f )]}-V 2 ^- (13) 
The graph of u versus H and xm versus H can be seen in fig. 2. We can select the 
modulation frequency u to control the fix points in periodic optical potential. Here we set 
dimensionless parameter u = 4.34, which corresponds to 3.67KHz. 

We use the symplectic integrators |T(| , [[□]] to solve the equations of motion because the 
Hamiltonian evolution preserves the Poisson bracket relation {x(t),p x (t)} = 1 

We plot the stroboscopic potrait of the system at times t = (2tt/uj)s, where s is an integer 
referred to as the strobe number. From figure 3,4,5 we can see that with the increase of e, 
the motion of atoms will appear chaos. But there are some stable regions. An initial wide 
spectral distribution of atoms will result in some atoms trapped in these stable regions. 

Laser cooling and trapping techniques have the ability to cool the atom to very low 
velocities and trap them in very small region in momentum. Therefore the appropriate de- 
scription of atomic dynamics is to use probability distribution on phase space (x,y,p x ,p y ). 



We define a classical state to be a probability measure on phase space of the form 
Q(x,y,p x ,Py)dxdydp x dp y , the density of probability satisfies the Liouville equation 

§ = {^W (14) 

where{, } is the Poisson bracket. The equation can be solved by the method of charac- 
teristics. To simulate the experiment, we assume initially atoms are umiformly distributed 
on \x\ < C and \y\ < C region, where C is a constant choosen to ensure the major fixed 
points are included. The momentum distributions for p x and p y are Gaussian distributions. 
Therefore 

Qo(x,y,p x ,p y ) = Qo(x)Qo(y)Qo(p x )Qo(p y ), (15) 

where 

QM = exp [-(p x - p x (0)) 2 /2a Pi }. (16) 

The variances of p x and p y are related to the temperature Tj 

a Pi = ksT./Pl (17) 



The two dimensional symplectic integrators |10| are used to keep the Poisson bracket 
relations during computation 

{q i ,p j } = S ij (18) 

Fig. 6 shows for optical potential when it is no modulation, atoms will accumulate 
around the fix pointx = y = (fig. 3). when the modulation is added the atoms will diffuse 
but will accumulate around several rings. With the increase of modulation amplitude, more 
atoms accumulate around rings and less atoms around x — y — (fig. 7,8). The variances 
a Px and a Py are taken as 0.05 in computation, which corresponds to temperature T x and T y 
approximately recoil temperature. If the variance is narrower rings become clearer. 
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IV. CONCLUSION AND DISCUSSION 



In summary, we have shown that it has chaotic dynamics for mudulated far-off-detuning 
donut beam. For atomic momenta p x , p y which have Gaussian distributions, atoms will be 
trapped in rings when the optical potential is modulated. If at some moment the optical 
potential is withdrawn, atoms will expand freely and will keep the shape of rings becasue 
the momentua is symmetric in (x,y). Therefore the two dimensional atomic distribution in 
(x, y) can be detected using TOF technique. 

Ackowlegement: Dr. S. Dyrting at Hong Kong University of Science and Technology 
gave one of authors helpful advice about two dimensional simplectic integrators. 
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FIGURES 




FIG. 1. The relative intensity of donut beam for charge 1 = 2, where w is the beam waist. 
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FIG. 2. The frequency of motion u>o ~ H and xm ~ H. 
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FIG. 5. Stroboscopic potrait of the system with e = 0.7, p x (0) = 0,p y = 1 and y = O.The 
maximun strobe number is 500. 
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FIG. 6. The atomic distribution in (x,y) plane at the strobe number 50 for e = 0. The 10000 
points were taken in phase space. The atoms were initially distributed on \x\ < 0.3 and \y\ < 0.3 
region. The momenta of p x , p y are Guassian distributions and cr Px = a p = 0.05. 
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FIG. 7. The atomic distribution in (x,y) plane at the strobe number 50 for e = 0.5. The 10000 
points were taken in phase space. The atoms were initially distributed on \x\ < 0.3 and \y\ < 0.3 
region. The momenta of p x , p y are Guassian distributions and o~ Px = o~ p = 0.05. 
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FIG. 8. The atomic distribution in (x,y) plane at the strobe number 50 for e = 0.7. The 10000 
points were taken in phase space. The atoms were initially distributed on \x\ < 0.3 and \y\ < 0.3 
region. The momenta of p x , p y are Guassian distributions and o~ Px = o~ p = 0.05. 
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